Frustration induced spin glass behaviour in Ba 3 RuTi 2 Og 
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Abstract 

Hexagonal Ba3RuTi20g consists of a layered triangular network of Ru 4+ (5 = 1) ions. Magnetic susceptibility x(T) and 
heat capacity data show no long range magnetic ordering down to 1.8 K. A Curie- Weiss (CW) fitting of x(T) yields a large 
antiferromagnetic CW temperature #cw = — 166 K. However, with low held, a splitting of zero held cooled (ZFC) and held 
cooled (FC) x(T) is observed below ~ 30K. Our measurements suggest that Ba3RuTi209 is a highly frustrated system and 
some parts of the material undergo a spin glass transition at 30 K. 

PACS numbers: 75.47.Lx, 75.50.Lk, 75.40.Cx 



I. INTRODUCTION 

In recent times, Ad and 5c? transition metal based ma- 
terials have been in the focus as their driving physics 
is different from their 3c? counterparts but have not 
been explored much. This variation arises because of 
their extended 'cf orbitals, a small onsite Coulomb en- 
ergy 'U', and a large spin-orbit coupling (SOC) [l|. 
The interplay between U and SOC makes these ma- 
terials interesting. Marry systems of the 4d/5d family 
like Sr 2 (lr/Ru/Rh)0 4 Ca9Ru0 4 0, Na 2 Ir0 3 @, 

Na 4 Ir 3 8 @, (Sr/Ca) 3 Ru 2 7 [all etc. have drawn at- 
tention in recent years. 

Recently, we reported results of our investigations of 
hexagonal Ba3ixTi20g with 5c? Ir 4+ (J e ff — 1/2) ions 
forming a 2D triangular network Q . A large site-disorder 
between Ir 4+ and Ti 4+ cations plays an important role 
in the magnetic properties of the sample. The system 
remains paramagnetic down to 0.35 K despite a large an- 
tiferromagnetic (AF) Curie- Weiss temperature 9qw ~ 
— 130 K. Magnetic heat capacity at low temperature fol- 
lows a power law with temperature. These results sug- 
gest that Ba3lrTi2 0g is in a spin liquid ground state with 
5c? Ir 4+ (Jeff = 1/2) in a 2D triangular network. The 
SOC presumably plays a significant role in determining 
the properties of the system. It will be interesting to 
study its 4d-analog (say, based on Ru) where spin-orbit 
coupling effect is expected to be of intermediate strength 
compared to that in 3d and 5c? elements. 

Dickson et al. [nj measured the magnetic susceptibil- 
ity variation with temperature of hexagonal Ba3RuTi20g 
in the limited temperature range 253 K to 323 K. Suscep- 
tibility data in a larger temperature range 77 — 335 K are 
reported in Ref. [llj. In both the papers a high, AF 6cw 

( 250 K) is inferred [H, OH]. These reports indicate 

that Ba3RuTi2 0g could also be a frustrated spin liquid 
candidate and a 4c? analog of Ba3lrTi2 0g with 5=1. 
Thus low temperature measurements on Ba3RuTi20g are 
warranted to confirm its magnetic ground state. 

In this paper we report preparation, structural anal- 
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ysis, magnetic susceptibility and heat capacity measure- 
ments on Ba3RuTi20g. The system crystallizes in the 
space group P63inc with a large site sharing between 
Ru 4+ and Ti 4+ ions as reported earlier [l2| and similar to 
that in Ba3irTi 2 0g [j|. Susceptibility data show no long 
range ordering down to 2K but splitting between zero 
field cooled (ZFC) and field cooled (FC) data is found 
below 30 K. A large 6> C w (AF) obtained from CW fitting 
indicates strong correlation between the magnetic ions. 
No anomaly, even at 30 K, is found in heat capacity mea- 
surements down to 1.8 K. These measurements suggest 
that the system is highly frustrated and perhaps a frac- 
tion of the spins participate in a spin glass state below 
30K. 



II. EXPERIMENTAL DETAILS 

A polycrystalline sample of Ba3RuTi2 0g was prepared 
by conventional solid state reaction method using high 
purity BaC03, Ti02 and Ru metal powder as starting 
materials. Stoichiometric amount of the starting mate- 
rials were mixed thoroughly, pressed into a pellet and 
calcined at 900°C for 15 h. After calcination, the pellet 
was crushed into powder, pelletized and fired at 1100°C 
for 100 h with several intermediate grindings. 

Powder x-ray diffraction (XRD) measurements were 
performed at room temperature with Cu K a radiation 
(A = 1.54182A) in a PANalytical X'Pert PRO diffrac- 
tometer using Si for calibration. Magnetization M mea- 
surements were carried out in the temperature T range 
2 - 400 K using Quantum Design PPMS and SQUID 
VSM. Heat capacity measurements were done in the tem- 
perature range 1.8 — 300 K and field H range — 9 T using 
the heat capacity attachment of Quantum Design PPMS. 



III. RESULTS AND DISCUSSIONS 

Single phase nature of our sample was confirmed from 
the XRD measurements. Maunders et al. have reported 
Ba3RuTi2 0g to crystallize in the space group P63TOC [13j. 
In this space group, the 2a site is occupied by the Ti 4+ 
10ns and in the 2b sites Ti 4+ /Ru 4+ ions reside in an or- 
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Figure 1: (i) Ideal structure of Ba3RuTi2C>9 without any site 
disorder between Ti 4+ and Ru 4+ ions. The Ru 4+ spins form 
an edge shared triangular network parallel to the ab plane as 
shown, (ii) RuTiOg bioctahedra are shown (iii) A possible 
arrang ement of Ti 4+ and Ru 4+ ions in the ab plane is shown 
when about 33% of Ru 4+ ions exchange site with Ti 4+ ions 
in the 2b site. 



dered manner. In the ideal case, the ordered arrange- 
ment of Ti 4+ and Ru 4+ ions form face-sharing RuTiOg 
bioctahedra (Fig. [TJii) ) and the Ru 4+ spins form an 
edge-shared triangular network in the ab plane (as shown 
in Fig. [TJi) ) . However, a site sharing between Ti 4+ 
and Ru 4+ ions is expected as the two ions are of sim- 
ilar ionic size. In Ref. four possible arrangements 
of Ru 4+ /Ti 4+ ions have been suggested in BaaRuT^Og. 
Among them two are most probable where RuTiOg and 
Ru20g/Ti20g bioctahedra are formed. We refined of our 
XRD data using the space group P63?tjc and considering 
a site sharing between Ti 4+ and Ru 4+ ions. We have 
found almost 31% site exchange of Ru 4+ ions with Ti 4+ 
ions in the 2b site and an 8% site sharing with 2a site 
Ti 4+ ions. This situation is very similar to that found in 
Ref. pjj]. In Ba3ixTi20g as well, we found a large site 
sharing of Ir 4+ ions with Ti 4+ ions both in 2b as well as 
2a sites |9j. 

Fig. [5] shows the XRD refinement pattern which yields 
R p = 5.7%, R wp = 7.7% and goodness of fit (GOF)= 
9.700. The lattice parameters obtained from the refine- 
ment are shown in Table U and compared with the val- 
ues with Ba3RuTi20g (reported earlier), Ba3ixTi20g and 
BaTi03 (parent compound). Values obtained by us are 
in good agreement with others. The atomic positions and 
occupancies resulting from our refinement are shown in 
Table IE 

As a result of site sharing between Ti 4+ and Ru 4+ ions, 
the Ru 4+ (5 = 1) edge shared triangular plane will be 
highly depleted. A probable situation is shown in Fig. 
[IJiii) where many of the misplaced Ru 4+ ions can still 
manage to interact with other Ru 4+ ions forming clusters, 
while some, called orphan spins, remain isolated. 

Further we have measured susceptibility of the mate- 



Table I: Lattice constants of Ba3RuTi 2 09 obtained from re- 
finement and comparison with others. 



Material 


Space group 


a (A) 


c(A) 


Source 


Ba 3 RuTi 2 9 


P6 3 mc 


5.7204 


14.0109 


this work 


Ba 3 RuTi 2 9 


P63771C 


5.7056 


14.0093 


Ref. [13] 


Ba 3 IrTi 2 9 


P6smc 


5.7214 


14.0721 


Ref. [9J 


BaTiOa 


P63 /mmc 


5.7238 


13.9649 


Ref. [14] 



Table II: Atomic parameters obtained by refining x-ray pow- 
der diffraction data for Ba3RuTi20g at room temperature 
with the space group P63mc. 





X 


y 


z 


g 


Ba(l) 


2a 








0.249 


1.00 


Ba(2) 


2b 


1/3 


2/3 


0.085 


1.00 


Ba(3) 


2b 


1/3 


2/3 


0.397 


1.00 


Ru(l) 


2b 


1/3 


2/3 


0.65 


0.61 


Ti(l) 


2b 


1/3 


2/3 


0.65 


0.39 


Ti(2) 


2b 


1/3 


2/3 


0.839 


0.69 


Ru(2) 


2b 


1/3 


2/3 


0.839 


0.31 


Ti(3) 


2a 








0.510 


0.92 


Ru(3) 


2a 








0.510 


0.08 


0(1) 


6c 


0.16098 


0.83898 


0.57538 


1.00 


0(2) 


6c 


0.48859 


0.51138 


0.74989 


1.00 


0(3) 


6c 


0.16098 


0.83898 


0.91548 


1.00 



rial at different fields as shown in Fig. [3J No long range 
magnetic ordering down to 2 K is found whereas a split- 
ting between ZFC and FC data at 30 K is seen which may 
be indicative of a transition to a glassy state. The split- 
ting is large, which indicates that a reasonable fraction 
of the spins are participating in the glassy phase. With 
increasing field, the splitting becomes less prominent and 
vanishes at an applied field of 10 kOe. 

We have fitted susceptibility data with CW formula 
X = Xo + C/ (T — 9) in the temperature range 100 — 350 K 
(shown in Fig 0]) which yields the temperature indepen- 
dent susceptibility xo = 2.05 x 10~ 4 cm 3 /mol, Curie con- 
stant C = 0.7cm 3 K/mol and #cw = — 166 K. For insu- 
lating oxides one can consider xo = Xcore + Xw, where 
Xcore is the core diamagnetic susceptibility and Xw is 
Van-Vleck paramagnetic susceptibility. For our sample 
Xcore — —2.32 x 10~ 4 cm 3 /mol [l5j|, which results in 
Xw = 4.37xl0~ 4 cm 3 /mol. This value of Xw is simi- 
lar to that found in other Ru based insulating oxides like 
La 2 Ru0 5 [3, La 2 LiRu0 6 [H etc. Inverse susceptibility 
is shown on the right axis of Fig. 0] (after subtracting the 
temperature independent part xo) which is linear with 
temperature in a wide range of temperature and devi- 
ates from linearity below ~ 80 K. 

We now compare our susceptibility data with that of 
Ref. [lH by Bryne et al. A slope change was found by 
Bryne et al. in the inverse susceptibility of Ba3RuTi20g 
at 150 K which is absent in our data. Moreover, they 
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Figure 2: Refinement pattern of Ba3RuTi2C>9 with space 
group P63771C. Light blue circles indicate experimental pat- 
tern and the red line indicate the theoretical pattern while 
the green line is the difference between the two. The Bragg 
positions are shown as vertical lines. The (hkl) planes corre- 
sponding to the peaks are also marked in the figure. 
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Figure 3: Magnetic susceptibility of Ba3RuTi2C>9 at different 
fields is shown as a function of temperature in semi-log scale. 
Solid symbols represent FC and open symbols represent ZFC 
measurements. 



found #cw = — 256K and fi e ff — 2.86/is which are 
higher than 8cw = — 166 K and fJ* e ff = 2.37/is obtained 
from our analysis. To clarify this discrepancy, we have 
reanalysed the published data as shown in Fig. 21 The 
green triangles represent the inverse susceptibility data 
as published in Ref. Corresponding susceptibility 

is shown as blue circles on the left axis and fitted by us 
with the CW law in the range 118 — 335 K. This fitting 
gives xob = 5.55 x 10~ 4 cm 3 /mol, C = 0.54 cm 3 K/mol 
(f/> e ff = 2.08 Hb) and 8qw — — 116 K. These values are 
somewhat closer to the values obtained by us. Subtract- 
ing this xoB; we have plotted (x — Xob) 1 which is linear 



Figure 4: Magnetic susceptibility of Ba3RuTi20g and its fit- 
ting with the Curie- Weiss law (red solid line) . Inverse suscep- 
tibility (after subtracting \o) is shown as orange diamonds. 
Inverse susceptibility data of Ba3RuTi209 obtained from Ref. 
[IH is shown as green triangles. The corresponding suscepti- 
bility is shown as blue circles with its fit with CW law (violet 
solid line). Inverse susceptibility after subtracting xob (ob- 
tained from fit) is also shown as pink squares. The dashed 
line is a guide to eye. All susceptibilities are plotted on the 
left axis while inverse susceptibilities correspond to the right 
axis. The solid symbols denote the data from this work and 
open symbols are used for data from Ref. (Tl| . 

in the whole temperature range and the kink at 150 K 
is also absent. Apparently, Bryne et al. did not use the 
correct xo leading them to infer an incorrect Curie term. 
Also, they calculated fi e ff based on a single suscepti- 
bility data point (one temperature) which certainly can 
mislead. 

The C value obtained from CW fit of our data is al- 
most 70% of the paramagnetic 5 = 1 value. Normally 
in 3d transition metal oxides the C value turns out as 
expected from a pure spin contribution. But for 5d tran- 
sition metal oxides the C value obtained is normally much 
smaller than expected from a pure spin contribution. In 
case of Ba3lrTi20g this was only 40% of the expected 
spin-only value [Jj- For Ad based systems the orbital 
effect is expected to be in between 3d and 5d materi- 
als. Hence the reduction in C obtained from our analysis 
is most likely an effect of spin-orbit coupling present in 
the system. A high (— 166 K) #cw obtained from CW 
fit confirms that the Ru 4+ spins are strongly correlated. 
Probably a large site sharing is responsible in forming 
magnetic clusters which give rise to the glassy behavior 
below 30K. 

Heat capacity (Cp) of a material gives prominent sig- 
natures pertaining to any structural or magnetic transi- 
tions. We have measured Cp at different fields as shown 
in Fig. [5] Surprisingly, no anomaly, even at 30 K is found 
in Cp measurements. In the low temperature region, Cp 
depends on applied field which can be seen clearly in 
the inset of Fig. [3] showing Cp /T as a function of T. 
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Figure 5: Heat capacity of Ba3RuTi20g measured at different 
fields is shown in a log-log scale. Inset: Cp/T is shown as a 
function of T in the low-T region. 



This field dependence could be arising from a Schottky 
anomaly of orphan spins present in the sample. The ab- 
sence of an anomaly at the spin glass transition tempera- 
ture (30 K) could be due to a large contribution from the 
crystal part supressing a small magnetic part at such high 
temperature. In this context we should mention that no 
anomaly in heat capacity is found in case of Ba3lrTi20g 
[| and Ba 3 TiIr 2 9 [lj| where ZFC and FC susceptibili- 
ties split up at 80 K and 60 K, respectively. Also, it is pos- 
sible that only a fraction of the spins are involved in the 
gjassy phase while the rest are not, similar to Ba3lrTi20g 

The total heat capacity of the sample has three contri- 
butions namely, (i) lattice contribution (C a t) (ii) mag- 
netic contribution (Cm) and (iii) Schottky contribution 
(Csch) from orphan spins. The nuclear Schottky anomaly 
normally appears at much lower temperature and prob- 
ably will not affect our data which is down to 1.8 K only. 
To extract the magnetic heat capacity we have to sub- 
tract out the Ci a t and Cs c h from the total heat capacity. 
We now attempt to estimate the magnetic contribution 
by subtracting the other two parts. 

The generalised heat capacity expression for two level 
Schottky anomaly is [l9l |20| , 



C Sch (A) = R°° 
9i 



A 



exp 



(*st) 



1 



yi 



exp 



A 



(1) 



where A is the level splitting, R is the universal gas 
constant, fcg is the Boltzman constant and go and g\ 
are the degeneracies of ground state and excited state, 
respectively. For spin 5 = 1 (Ru 4+ ) systems, Eq. Q]cair 
be modified with the approximation as go = 1 and gi = 2 

We have subtracted C P (0T) from C P (9T) data and 
further divided that by temperature. Any lattice con- 



Figure 6: Magnetic heat capacity of Ba3RuTi20g at 0T and 
9T are shown in a log-log scale. The solid lines shows power 
law with temperature with the power shown in figure. Inset: 
[C P (9T)- C P (0T)]/T and its fit with Schottky model (solid 
line) is shown. 



tribution is therefore removed by this procedure. The 
[Cp(9T)-C P (0T)]/T data are shown in the inset of Fig. 
E] with its fit to /[C S ch(A 9T ) - C Sch (A 0T )}/T . In this 
expression, / is the percentage of orphan spins present in 
the sample and Csch(AgT), and Csch(AoT) are the Schot- 
tky anomalies corresponding to the level splittings of AgT 
and Aot with applied magnetic fields 9T and 0T, respec- 
tively. The fit is reasonably good and suggests that the 
field dependence of Cp is coming from Cs c h only. From 
the fit we found A 9T /k B = 27. 7K, A 0T /k B = 1.5 K and 
~ 1% orphan spins are present in the system. Putting 
the corresponding A (obtained from fitting) in Eq. Q]and 
multiplying by /, we can get the Schottky contribution 
(Csch) to the heat capacity for a particular field. 

To subtract the lattice contribution (Ci a t) from the 
total heat capacity (Cp) one needs a non- magnetic ana- 
log of the system. Unfortunately, we didn't find a good 
non- magnetic analog of Ba3RuTi2 0g. However we have 
earlier used Ba 3 ZnSb2 0g (taken from Ref. [2l|) as 
a non-magnetic analog in the case of Ba 3 IrTi 2 09 @. 
Ba3ZnSb20g can be used as a non magnetic analog in 
the present case as well since both have similar crystal 
structure. The heat capacities of these two samples have 
a mismatch at higher temperature due to the difference in 
their molecular mass and volume. We have used a multi- 
fication factor of ~ 0.67 so that Cp data of Ba3ZnSb20g 
matches with that of Ba3RuTi2 0g in the nominal tem- 
perature range 40 — 60K. The corrected heat capacity 
of Ba3ZnSb20g can be considered as the lattice contri- 
bution (Cat) to the total heat capacity of Ba3RuTi20g. 
This Ci a t is subtracted from Cp and we are left with the 
Cm and Cs c h- Further we have subtracted Cs c h obtained 
from our previous analysis and finally obtained Cm- 

The Cm is independent of field as shown in Fig. [5] 
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and follows a power law in temperature with power 0.6 
in the range 4.3 — 10 K and with power 2.1 below 2.7 K. 
This suggests that spin excitations are present in this 
system. We have also estimated the entropy change 
ASm=0.36 J/mol-K by integrating Cm/T. This value is 
much smaller then expected value of 9.13 J/mol-K. This 
reduction could be partly due to a strong spin-orbit cou- 
pling present in the system. Most likely other reasons 
like error in estimation of Cut due to lack of a proper 
non-magnetic analog are also responsible. 

IV. CONCLUSIONS 

We have prepared Ba3RuTi20g (a 4d analog of the 
5d system Ba3irTi2 0g) and measured its magnetic sus- 
ceptibility and heat capacity at different fields. No long 
range magnetic ordering is found down to 1.8 K despite 
a large # cw obtained from CW fitting of susceptibility 
data. This indicates that the spins are strongly correlated 



and the system is highly frustrated. Magnetic moments 
are reduced probably because of a spin-orbit coupling 
present in this system. Splitting of ZFC and FC sus- 
ceptibility data at 30 K indicates formation of a glassy 
phase. However, absence of an heat capacity anomaly 
at this temperature probably suggests that a fraction 
of the spins participate in this glassy phase. Magnetic 
heat capacity of this material is independent of field and 
approaches a power law at low temperature with power 
2.1. These measurements suggest that Ba 3 RuTi 2 09 is 
a highly frustrated system but the nature of its ground 
state is not very clear. Local probe investigations like 
^iSR could be useful to check if the spin freezing happens 
in the bulk of the sample or not. 
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